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Abstract 

If the observation by OPERA of apparently superluminal neutrinos is correct, the 
Lagrangian for second-generation leptons must break Lorentz invariance. We calculate 
the effects of an energy-independent change in the neutrino speed on another observable, 



the charged pion decay rate. The rate decreases by a factor 

where (vu) is the (directionally averaged) neutrino speed in the pion's rest frame. This 
provides a completely independent experimental observable that is sensitive to the same 
forms of Lorentz violation as a neutrino time of flight measurement. 
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The recent announcement by the OPERA experiment that muon neutrinos may be 
travehng faster than hght with speeds of ~ 1 + (2.5 x 10^^) [Ij has attracted a tremen- 
dous amount of attention. Whether this experimental conclusion will withstand further 
scrutiny remains to be seen. However, the mere possibility that superluminal travel might 
occur has already produced an explosion of interest in theories that do not precisely obey 
Lorentz symmetry. 

Before the intriguing neutrino results were reported, there was already a thriving com- 
munity of physicists working on such Lorentz- violating theories. This community includes 
both experimenters and theorists. One of the most important theoretical developments 
in the field was an effective quantum field theory capable of describing all possible forms 
of Lorentz violation that might occur at low energies with standard model fields [21 [3] . 
This theory is known as the standard model extension (SME). Both the renormalizabil- 
ity [H El E] and stability of the SME have been studied. The SME is now the standard 
framework for parameterizing the results of experimental Lorentz tests. Many different 
systems and phenomena may be used for such tests, and up-to-date information about 
most of the resulting constraints may be found in [8]. 

In this paper, we shall use the SME formalism to examine one very important conse- 
quence that the existence of superluminal neutrinos could have. If the neutrino energy- 
momentum relation is modified, the rates of scattering and decay processes involving these 
neutrinos will necessarily be modified. Of particular interest are meson decays. The phase 
space available to the decay products in a process such as vr"*" — )■ /x"*" -|- depends crucially 
on the neutrino velocity. This provides a completely new observable, which is nonetheless 
sensitive to some of the same parameters as the OPERA experiment. Precision decay 
measurements can thus be directly compared with the results of a neutrino time of flight 
experiment. 

We shall calculate how an energy-independent change in the muon neutrino speed 
affects the decay rate of a charged pion. The choice of the pion is partially just for 
definiteness. However, we shall also find that the Lorentz-violating modifications to the 
decay rate are most pronounced for the lightest parent particles. While this makes pions 
the ideal species to study, our calculations naturally apply to any particle that decays via 
a W boson primarily to a muon and its associated neutrino. 

Moreover, we expect that similar phenomena would occur in more complicated decay 
processes, such as /i"*" — > e"^ + z/e + z/^. In fact, calculations of Cerenkov-like decay processes 
such as — )■ Ufj, + + e~ have already raised serious questions about the OPERA 
result. Such processes, which could be allowed for superluminal neutrinos, would have 
major effects on the neutrino beam in transit — effects which have seemingly not been 
observed PUni. 

In our calculations, we shall assume that the only Lorentz violation is in the lepton 
sector. The structures and interactions of the pion itself and the W that mediates the 
decay are exactly as in the standard model. This means that the entire calculation may 
be performed, as usual, in the rest frame of the pion. The kinematic factors associated 
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with the motion of the parent particle are unchanged from the conventional case. 

However, it is not possible in the SME formalism to have Lorentz violation for muon 
neutrinos but simultaneously to have a muon sector that is free of Lorentz violation. The 
lepton Lagrange density that contains a spin- and energy-independent modification to the 
neutrino velocity is 

L = iL{Y + c''^'-f^)D^L + iR{Y + c''^'^^)D^R + mass terms, (1) 

where L and R are the left- and right-chiral lepton multiplets, 

R = [in] . (2) 

and Dn is the covariant derivative, containing the electromagnetic and weak gauge in- 
teractions. SU{2)l gauge invariance requires that the Lorentz violation coefficients c'"'^ 
be the same for the charged lepton and its neutrino. (Actually, it is possible to have a 
different set of coefficients for the right-chiral charged lepton field, but those coefficients 
are irrelevant in reactions involving charged current weak interactions.) The c coefficients 
together form a traceless background tensor, which imbues spacetime with a preferred 
frame structure. While the c coefficients control violations both of isotropy and of boost 
invariance, we shall mostly be concerned here with the isotropic element of the tensor cqq. 
Because Lorentz violation is a small effect, we may neglect all terms beyond first order in 
c. 

The choice of an energy-independent modification of the lepton velocities is motivated 
by two important factors. First, the OPERA data does not show evidence of any energy 
dependence in the neutrino speed. Second, the observed effect is rather large, suggesting 
that the modification of the velocity is not an effect suppressed by any powers of a large 
energy scale. This requires that the Lorentz-violating operator responsible for the effect 
be renormalizable, leading to the choice of c. 

There are very few direct constraints on the c coefficients for second-generation leptons. 
The absence of the photon decay process 7 — )■ fi~^ + fi~ up to TeV photon energies 
constrains the possibility that the leptons might have limiting velocities less than the 
speed of light pjj, but there are few constraints on the possibility that muons and their 
corresponding neutrinos might move faster than light. Superluminal muons would emit 
vacuum Cerenkov radiation [12], but to constrain the muon c coefficients at the level seen 
by OPERA would require careful observation of ~ 15 GeV muons. 

The strongest bounds on other forms of muon Lorentz violation are based on (7^ — 2 
measurements This data could also be used to constrain the muon c, potentially at 
the 10~^ level. The behavior of a fermion in an external field, with a form of photon-sector 
Lorentz violation that is essentially equivalent to c, has been calculated. The result is 
that the spin couples to both the electric and magnetic fields in a potentially anisotropic 
fashion [14j. A full analysis of this phenomenon would need to include the effects of the 
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substantial E fields present in muon — 2 experiments as well as the possible influences 
of other SME parameters. 

On the other hand, there are some constraints on c-type Lorentz violation in the 
neutrino sector. If the c coefficients are nondiagonal in flavor space, they will lead to neu- 
trino oscillations, with a different characteristic energy dependence from the oscillations 
generated by nondiagonal neutrino masses [15]. This possibility has been constrained, 
and Lorentz violation as the principal cause of observed neutrino oscillations has already 
been ruled out [16]. Oscillation measurements can also constrain the differences between 
the flavor- diagonal c parameters for different species, since if the different neutrino flavors 
had substantially different speeds, the wave packet overlap necessary for oscillations could 
be destroyed. With only isotropic c coefficients, the oscillation data require a very high 
degree of flavor independence (at the 10~^^ level), which would make the OPERA data 
incompatible with numerous measurements of electron-sector Lorentz violation [T71 [18] . 
However, the number of SME parameters that can influence neutrino oscillations is very 
large, and it is not clear whether this conclusion applies in more general models. Finally, 
there is the famous time of flight measurement made on neutrinos from SN 1987A [T^EU], 
which (assuming the interpretation of the data was correct) constrains the c coefficients 
for electron neutrinos at the 10~^ level. 

Constraints on the SME parameters are conventionally expressed in a sun-centered 
reference frame, with coordinates (X, Y, Z, T), with the Z-direction parallel to the Earth's 
rotation axis ^T\. The OPERA neutrino timing data was collected over a long period of 
time, effectively averaging over the range of neutrino directions v that are possible as the 
Earth rotates. As shown below in f|TT]) . the speed of a massless particle in the presence of 
the c parameters is 1 — cqo — C(^oj)Vj — CjkVjVk, where C(oj) is the symmetrized combination 
%■ + ^jo- This makes the average speed of the OPERA neutrinos 1 — Cqo — C(oz) cos 9 — 
\{cxx + cyy) sin^ 6 — czz cos^ 6*, where 6 is the angle between the CERN-to-Gran Sasso 
direction and the Z-axis, for which cos^^ ~ —0.41. 

We now turn to the calculation of the pion decay rate. This rate is not directly sensitive 
to all the c parameters that might have affected the average neutrino speed. Because, 
in the absence of Lorentz violation, the distribution of the decay products is isotropic, 
the modifications to the total decay rate can depend only on the isotropic Cqo = Cjj, as 
measured in the decaying particle's rest frame. We may therefore treat c as if it were a 
diagonal tensor, with elements [cqo, |coo, |coo, |coo] in the center of mass frame. If this 
isotropic c is the sole contributor to the neutrinos' faster-than-light travel, v^^ = 1 — |coo. 
However, the value of Cqo will depend on the pion momentum, and other components of 
c will affect the differences in the decay rates of pions with different velocities. 

Although the ultimate form taken by the modified decay rate is quite simple, the 
calculations leading up to it are somewhat tricky. Methods for calculating cross sections 
and decay rates in the SME context have been worked out [22]. The formulas for physical 
observables can be put in forms quite similar to those that appear conventionally. How- 
ever, in these formulas, both the matrix element for the decay and the kinematics are 
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modified from tfieir conventional forms. Tfie kinematic modifications, wliicli are related 
to clianges in tlie pliase space available to the external particles, are often more important 
than the changes to the matrix element. 

The matrix element Ai can be determined from the same tree-level Feynman diagram 
that governs conventional pion decay. Since there is assumed to be no Lorentz violation 
in the meson or W boson sectors, only those parts of the diagram involving the final 
muon and neutrino need to be considered. This means that A4 may be modified in two 
ways. The first is via the replacement of 7^ at the muon-neutrino-gauge boson vertex 
with 'y'^ + c'^^^u. The free index fi is ultimately to be contracted with the total momentum 
vector p^, which has only a time component. This amounts to a rescaling of the matrix 
element simply by a factor of (1 + cqo)- 

The second modification of the matrix element comes from the modified spinors, which 
obey the modified Dirac equation 



1 + coo)7o^ 



1 



Coo IjPj 



— m 



u{p) = 0. 



(3) 



The spinor solutions are standard, except for an evident rescaling of E and pj. Therefore, 
if they are normalized to obey u** {p)u''{p) = ^ESss' [and analogously for v{p)], the closure 
relations needed for calculating cross sections are u'^{p)u'^{p) = 7°-E+ (l — |coo) 7jPj + 
m and v'^{p)v^{p) = 7°-E'+ (l — |coo) 'JjPj —m. This makes the matrix element squared 
proportional to 



\M\^ oc (l + Coo)'5^tr |t;^'(-PM)7o(l-75)<(p;)«M(^M)(l + 75)7o<(-p;)} (4) 



;i + Coo)'tr 



3 Coo ) lj{-P,.)j 



7o(l - 75) 



3 Coo ) lj{Pt^)J 



;i + 75)70 



2(1 + -Coo ) Pti 



E„ 



1 - 3C00 ) 



(5) 
(6) 



In these formulas, and Pi, denote the magnitudes of the muon and neutrino three- 
momenta. (Despite the indices, these symbols do not represent four- vectors.) E^ and E^, 
are the corresponding energies. The derivation of (EI) used the facts that the neutrino 
momentum p^, is opposite to the muon momentum p)j and that E^ = (^1 — |coo) Pu, which 
is an obvious consequence of the massless Dirac equation. 

Further simplification requires analysis of the more complicated dispersion relation for 
the muon. The energy of a fermion, including the leading order effects of c, is 



E 



m? + PjPj){l - 2coo) - ^CjkPjPk - 2cojPj^/m^ + pjPj, 



(7) 
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which reduces to E = y {1 — 2cQo)m? + (1 — ^Coo)pjPj when only the isotropic cqo is con- 
sidered. Energy and momentum conservation then dictate that the muon momentum is 



2(2+g^) 
3(l-€2) 



Coo 



, where ^ = — is the muon-pion mass ratio — the main dimen- 



sionless parameter determining the kinematics of the Lorentz-invariant case. The energies 



2e 



Coo ) and E„ 



1 - 



2(- 



yCoo 



of the decay products are E^, — —^m.,^ 

The quantity — (l — |coo) P,x is then ^^m^(l — 2coo), so the ratio of the matrix element 
squared in the presence of cqq to its value at cqq = is 



M 



Mr 



:Coo 



(8) 



with the subscript referring to the Lorentz-invariant values of particular quantities. 

As mentioned, the Lorentz violation changes not just the matrix element for this 
process, but also the phase space factors in the decay rate formula, since Cqo affects the 
free propagation of the daughter particles. The phase space integral appearing in the 
decay rate is 



(27r)3 2E^ (27r)3 2E, 

d^p,. 



167r2 



E^{p^)E^{p^) 



(27:f5\p,-p,-p,) 
6{m^ - E^- Ey). 



(9) 
(10) 



After eliminating p^,, it remains to evaluate the integral over the muon momentum, d^p^ — 

pfidpfidflfj,. The angular integration can be done by symmetry, yielding 47r; then the 

(5-function eliminates the dp^ integration, fixing at its physical value. The energy 5- 

-1 

function also introduces a factor -^{t^-k — — Ei) — \vij,-\-Vi,\'''^ ^ where the velocities 
are evaluated at the fixed momentum p^. 

The neutrino velocity is simple, f = 1 — |coo, 
general expression for a massive particle's group velocity is 



but the muon case is trickier. The 



Cqg- 



which, in the isotropic limit, reduces to 



2c,- 



+ 



+ Ckl- 



TTjTrkTTl 



+ 7r2)3/2 



2c, 



Oj, 



(11) 



Pi 



3 3(m2+p2)3/2 



Coo 



:i2) 



Evaluated at the real muon momentum, this becomes 



1-e 



3(i+e)Vi-e i+e 



5 
3 



Coo 



(13) 



5 



which makes + Vy the considerably simpler 

2 



3(1 +e 



-Coo 



(14) 



It is now possible to assemble all the pieces of the cross section. Since there is no 
Lorentz violation in the meson sector, the kinematic factors associated with the initial 
state are not modified. The remaining effects give 



r 



1 + 



E 



E„ 



M 



Mr 



1-C 



(15) 
(16) 



The first factor on the right-hand- side of f|T5l) arises from the in d^Pfj,. 

The final formula f|T6|) is quite simple, and it depends on the same single quantity 
1 — that typically characterizes leptonic decay rates. [The simplicity of the final result 
also suggests that there may be a more convenient way of organizing this calculation, in 
which complicated rational functions of such as those in f[T^ . need not appear.] Using 
the physical value of ^ = 0.757 for pion decay, T = ro(l + 9.4coo)- For heavier mesons, the 
dependence on coo will be smaller. Note that coo can also be expressed as —j{{vu) — 1), 
where (viy) is the neutrino speed (averaged over all directions v) in the rest frame of the 
decaying pion. 

The modification to F given in ( !T6l) depends on the mass of the parent particle. Most 
of the c dependence in F comes from the kinematic factors rather than from the dynamical 
matrix element. The fact that the kinematics are more strongly affected when ^ is close 
to 1 is relatively straightforward to understand. At small momenta, a small addition to 
the velocity produces a much larger fractional change in the available phase space than 
would a similar term at large momenta. This tends to make the decay rate more sensitive 
to c when the energy available to the decay products is small. 

Since superluminal neutrinos correspond to negative values of Coo, superluminal behav- 
ior will tend to reduce the rate at which a pion will decay. This is natural, since the decay 
products (both the neutrinos and the corresponding muons) have dispersion relations for 
which the energy increases more rapidly than usual as a function of the three-momentum. 
Since the energy available in the decay is fixed to be m.„, the negative value of cqo decreases 
the momentum that the daughter particles can carry and consequently also decreases the 
available phase space. 

There are two straightforward kinds of comparisons that might be made, in order to 
test the hypothesis of a nonzero muon neutrino c. The first involves a comparison of the 
decay rates for pions (or other mesons) moving with different velocities. The Coo in ( !T6|) 
is the coefficient in the decaying particle's rest frame. In terms of the coefficients ctt and 
ctj in the sun-centered frame, the Coo is 



Coo = ll [cTT + C(TJ){v-„)j + CjK{V-n)j{Vn)K\ 



(17) 
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where is the meson velocity in the sun-centered frame, and is the corresponding 
Lorentz factor 7 = (1 — vl)'^^"^. 

Generally, the size of the Lorentz violation effect grows as the square of the meson 
energy. (This novel phenomenon is, of course, in addition to the conventional time dilation 
effect that lengthens the apparent lifetime of the meson, as measured by a stationary 
observer, by Although the absolute magnitude of the pion decay rate depends on 
difficult-to-compute hadronic effects, comparisons of the decay rates for mesons with 
different boosts can be used to search for a nonzero ctt- Most obviously, measurements 
of highly boosted decays may be compared with the reference value of the decay lifetime for 
a meson at rest. Moreover, studying the decays of mesons moving in different directions 
adds potential sensitivity to the other coefficients C(tj) and cjk- Even if, at a given 
source facility, the mesons produced are traveling in relatively tight beam, oriented along 
a fixed compass direction, the rotation of the Earth will make it possible to sample 
multiple directions v. The current best value of the charged pion lifetime is accurate at 
the 2 X 10~^ level in measurements with nonrelativistic pions. Experiments of comparable 
precision with relativistic pions with ^ 3 would be sensitive to anomalies in the neutrino 
dispersion relation at the level reported by OPERA. At higher energies, experiments can 
be less precise by a factor of 7^ and still be capable of confirming the novel effects. 

The other way in which this effect might be straightforwardly tested involves com- 
parisons of multiple meson decay modes. Light charged mesons are generally much more 
likely to decay into muon and neutrino pairs than electron and neutrino pairs. However, 
the latter decay does occur, and the ratio of the two decay rates can be predicted quite 
precisely in the standard model. Tests of lepton universality, which measure the relative 
sizes of the branching ratios for decays such as vr^ — ?■ /x"^ + z/^ and vr"'' — e"*" + z/g, are 
therefore sensitive to the muon c coefficients. 

Of course, the tt"*" — )■ e"*" + v^. decay rate is sensitive to the c coefficients in the electron 
sector. However, the coefficients for electrons (and, therefore, their associated neutrinos) 
are much more tightly constrained than the coefficients for the second-generation leptons. 
The reason is that electrons are stable and extremely numerous. They play a crucial role 
in many high-energy phenomena, and their energy-momentum relation has been mapped 
out very carefully, even up to PeV energies. This can be done by studying the radiation 
emitted by highly boosted electrons; such electromagnetic emissions are sensitive probes 
of both electron velocities and energies. By looking at the ways various processes would 
be modified by the inclusion of electron c coefficients and comparing these predictions 
to observations of energetic electrons — both in astrophysical sources [23] and at particle 
accelerators [21] — all the relevant coefficients may be constrained at the 10~^^ level or 
better. 

Finally, it is worth noting that for sufficiently large boosts, the decay vr"*" — )■ /i"*"-!-//^ may 
be impossible. If > ^ {ml - m2)/2[-CTT - qTj){v^)j - cjk{v^)j{v^)k], the greater- 
than-normal growth of the muon and neutrino energies makes the decay energetically 
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impossible [H |25]. The pion energy 77rm^ is insufficient to produce the highly boosted 
decay products. In contrast, in a purely phenomenalistic model, which violates gauge 
invariance by giving only the neutrino a modified dispersion relation, there is no threshold 
above which the pion is stable [2B]. 

If neutrinos do move faster than light, this will have many consequences in particle 
physics. Some of these consequences are rather counter-intuitive. For example, without 
also violating electroweak gauge invariance, it is impossible to have an energy-independent 
form of Lorentz violation for neutrinos without having it for the charged leptons as well. 
A complete dynamical calculation of the pion decay rate shows that it is fairly sensitive 
to this kind of Lorentz violation, and this provides a completely independent avenue for 
testing the fascinating OPERA result. 
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